Search for: All records

Mass segregation is seen in many star clusters, but whether massive stars form in the center of a cluster or migrate there dynamically is still debated.N-body simulations show that early dynamical mass segregation is possible when sub-clusters merge to form a dense core with a small crossing time. However, the effect of gas dynamics on both the formation and dynamics of the stars could inhibit the formation of the dense core. We aim to study the dynamical mass segregation of star cluster models that include gas dynamics and selfconsistently form stars from the dense substructure in the gas. Our models use the TORCH framework, which is based on AMUSE and includes stellar and magnetized gas dynamics, as well as stellar evolution and feedback from radiation, stellar winds, and supernovae. Our models consist of three star clusters forming from initial turbulent spherical clouds of mass 10⁴, 10⁵, 10⁶M_⊙and radius 11.7 pc that have final stellar masses of 3.6 × 10³M_⊙, 6.5 × 10⁴M_⊙, and 8.9 × 10⁵M_⊙, respectively. There is no primordial mass segregation in the model by construction. All three clusters become dynamically mass segregated at early times via collapse confirming that this mechanism occurs within sub-clusters forming directly out of the dense substructure in the gas. The dynamics of the embedded gas and stellar feedback do not inhibit the collapse of the cluster. We find that each model cluster becomes mass segregated within 2 Myr of the onset of star formation, reaching the levels observed in young clusters in the Milky Way. However, we note that the exact values are highly time-variable during these early phases of evolution. Massive stars that segregate to the center during core collapse are likely to be dynamically ejected, a process that can decrease the overall level of mass segregation again. 

Abstract We investigate how stellar feedback from the first stars (Population III) distributes metals through the interstellar and intergalactic medium using the star-by-star cosmological hydrodynamics simulation, Aeos. We find that energy injected from the supernovae (SNe) of the first stars is enough to expel a majority of gas and injected metals beyond the virial radius of halos with massM_dm ≲ 10⁷M_⊙, regardless of the number of SNe. This prevents self-enrichment and results in a nonmonotonic increase in metallicity at early times. Most minihalos (M_dm ≳ 10⁵M_⊙) do not retain significant fractions of the yields produced within their virial radii until they have grown to halo masses ofM_dm ≳ 10⁷M_⊙. The loss of metals to regions well beyond the virial radius delays the onset of enriched star formation and extends the period that Population III star formation can persist. We also explore the contributions of different nucleosynthetic channels to 10 individual elements. On the timescale of the simulation (lowest redshiftz= 14.3), enrichment is dominated by core-collapse supernovae for all elements, but with a significant contribution from asymptotic giant branch winds to thes-process elements, which are normally thought to only be important at late times. In this work, we establish important mechanisms for early chemical enrichment, which allows us to apply Aeosin later epochs to trace the evolution of enrichment during the complete transition from Population III to Population II stars. 

Abstract The Aeosproject introduces a series of high-resolution cosmological simulations that model star-by-star chemical enrichment and galaxy formation in the early Universe, achieving 1 pc resolution. These simulations capture the complexities of galaxy evolution within the first ~300 Myr by modeling individual stars and their feedback processes. By incorporating chemical yields from individual stars, Aeosgenerates galaxies with diverse stellar chemical abundances, linking them to hierarchical galaxy formation and early nucleosynthetic events. These simulations underscore the importance of chemical abundance patterns in ancient stars as vital probes of early nucleosynthesis, star formation histories, and galaxy formation. We examine the metallicity floors of various elements resulting from Population III enrichment, providing best-fit values for eight different metals (e.g., [O/H] = −4.0) to guide simulations without Population III models. Additionally, we identify galaxies that begin star formation with Population II after external enrichment and investigate the frequency of carbon-enhanced metal-poor stars at varying metallicities. The Aeossimulations offer detailed insights into the relationship between star formation, feedback, and chemical enrichment. Future work will extend these simulations to later epochs to interpret the diverse stellar populations of the Milky Way and its satellites. 

Abstract The first generations of stars left their chemical fingerprints on metal-poor stars in the Milky Way and its surrounding dwarf galaxies. While instantaneous and homogeneous enrichment implies that groups of conatal stars should have the same element abundances, small amplitudes of abundance scatter are seen at fixed [Fe/H]. Measurements of intrinsic abundance scatter have been made with small high-resolution spectroscopic data sets where measurement uncertainty is small compared to this scatter. In this work, we present a method to use mid-resolution survey data, which have larger errors, to make this measurement. Using APOGEE Data Release 17, we calculate the intrinsic scatter of Al, O, Mg, Si, Ti, Ni, and Mn relative to Fe for 333 metal-poor stars across six classical dwarf galaxies around the Milky Way, and 1604 stars across 19 globular clusters (GCs). We calibrate the reported abundance errors in bins of signal-to-noise ratio and [Fe/H] using a high-fidelity halo data set. Applying these calibrated errors to the APOGEE data, we find small amplitudes of average intrinsic abundance scatter in dwarf galaxies ranging from 0.03 to 0.09 dex, with a median value of 0.047 dex. For the GCs, we find intrinsic scatters ranging from 0.01 to 0.11 dex, with particularly high scatter for Al and O. Our measurements of intrinsic abundance scatter place important upper bounds, which are limited by our calibration, on the intrinsic scatter in these systems, as well as constraints on their underlying star formation history and mixing that we can look to simulations to interpret. 

Two main mechanisms have classically been proposed for the formation of runaway stars. In the binary supernova scenario (BSS), a massive star in a binary explodes as a supernova, ejecting its companion. In the dynamical ejection scenario, a star is ejected during a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway stars: the subcluster ejection scenario (SCES), where a subset of stars from an infalling subcluster is ejected out of the cluster via a tidal interaction with the contracting gravitational potential of the assembling cluster. We demonstrate the SCES in a star-by-star simulation of the formation of a young massive cluster from a 10⁶M_⊙gas cloud using theTORCHframework. This star cluster forms hierarchically through a sequence of subcluster mergers determined by the initial turbulent, spherical conditions of the gas. We find that these mergers drive the formation of runaway stars in our model. Late-forming subclusters fall into the central potential, where they are tidally disrupted, forming tidal tails of runaway stars that are distributed highly anisotropically. Runaways formed in the same SCES have similar ages, velocities, and ejection directions. Surveying observations, we identify several SCES candidate groups with anisotropic ejection directions. The SCES is capable of producing runaway binaries: two wide dynamical binaries in infalling subclusters were tightened through ejection. This allows for another velocity kick via subsequent via a subsequent BSS ejection. An SCES-BSS ejection is a possible avenue for the creation of hypervelocity stars unbound to the Galaxy. The SCES occurs when subcluster formation is resolved. We expect nonspherical initial gas distributions to increase the number of calculated runaway stars, bringing it closer to observed values. The observation of groups of runaway stars formed via the SCES can thus reveal the assembly history of their natal clusters. 

The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the TORCHframework, which uses the Astrophysical MUltipurpose Software Environment (AMUSE) to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgraded TORCHby implementing the N-body code PETAR, thereby enabling TORCHto handle massive clusters forming from 10⁶M_⊙clouds with ≥10⁵individual stars. We present results from TORCHsimulations of star clusters forming from 10⁴,  10⁵, and 10⁶M_⊙turbulent spherical gas clouds (named M4, M5, M6) of radiusR= 11.7 pc. We find that star formation is highly efficient and becomes more so at a higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities 2.325 × 10^1,2,3M_⊙pc⁻², after a free-fall time oft_ff= 6.7,2.1,0.67 Myr, we find that ∼30%, 40%, and 60% of the cloud mass has formed into stars, respectively. The end of simulation-integrated star formation efficiencies for M4, M5, and M6 areϵ_⋆ = M_⋆/M_cloud = 36%, 65%, and 85%. Observations of nearby clusters similar in mass and size to M4 have instantaneous star formation efficiencies ofϵ_inst ≤ 30%, which is slightly lower than the integrated star formation efficiency of M4. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense (Σ ≥ 10²M_⊙pc⁻²) and massive (M_cloud≥ 10⁵M_⊙). In such environments, stellar feedback from winds and radiation is not strong enough to counteract the gravity from gas and stars until a majority of the gas has formed into stars. 

Context.The star cluster initial mass function is observed to have an inverse power law exponent around 2, yet there is no consensus on what determines this distribution, and why some variation is observed in different galaxies. Furthermore, the cluster formation efficiency (CFE) covers a range of values, particularly when considering different environments. These clusters are often used to empirically constrain star formation and as fundamental units for stellar feedback models. Detailed galaxy models must therefore accurately capture the basic properties of observed clusters to be considered predictive. Aims.We study how feedback mechanisms acting on different timescales and with different energy budgets affect the star cluster mass function and CFE. Methods.We use hydrodynamical simulations of a dwarf galaxy as a laboratory to study star cluster formation. We test different combinations of stellar feedback mechanisms, including stellar winds, ionizing radiation, and supernovae (SNe). Results.Each feedback mechanism affects the CFE and cluster mass function. Increasing the feedback budget by combining the different types of feedback decreases the CFE by reducing the number of massive clusters. Ionizing radiation is found to be especially influential. This effect depends on the timing of feedback initiation, as shown by comparing early and late feedback. Early feedback occurs from ionizing radiation and stellar winds with onset immediately after a massive star is formed. Late feedback occurs when energy injection only starts after the main-sequence lifetime of the most massive SN progenitor, a timing that is further influenced by the choice of the most massive SN progenitor. Late feedback alone results in a broad, flat mass function, approaching a log-normal shape in the complete absence of feedback. Early feedback, on the other hand, produces a power-law cluster mass function with lower CFE, albeit with a steeper slope than that usually observed. 

We report the first high-resolution, detailed abundances of 21 elements for giants in the Galactic bulge/bar within 1° of the Galactic plane, where high extinction has rendered such studies challenging. Our high-signal-to-noise-ratio and high-resolution, near-infrared spectra of seven M giants in the inner bulge, located at (l,b) = (0°, +1°), are observed using the IGRINS spectrograph. We report the first multichemical study of the inner Galactic bulge by investigating, relative to a robust new solar neighborhood sample, the abundance trends of 21 elements, including the relatively difficult to study heavy elements. The elements studied are: F, Mg, Si, S, Ca, Na, Al, K, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, Ce, Nd, and Yb. We investigate bulge membership of all seven stars using distances and orbital simulations, and we find that the most metal-poor star may be a halo interloper. Our investigation also shows that the inner bulge as close as 1° north of the Galactic Center displays a similarity to the inner disk sequence, following the high-[α/Fe] envelope of the solar vicinity metal-rich population, though no firm conclusions for a different enrichment history are evident from this sample. We find a small fraction of metal-poor stars ([Fe/H] > −0.5), but most of our stars are mainly of supersolar metallicity. Fluorine is found to be enhanced at high metallicity compared to the solar neighborhood, but confirmation with a larger sample is required. We will apply this approach to explore the populations of the nuclear stellar disk and the nuclear star cluster.